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ABSTRACT 

We present newly calculated low-temperature opacities for gas with a primordial chemi- 
cal composition. In contrast to earlier calculations which took a pure metal-free Hydro- 
gen/Helium mixture, we take into account the small fractions of Deuterium and Lithium as 
resulting from Standard Big Bang Nucleosynthesis. Our opacity tables cover the density range 
-16 < logp [g cirT 3 ] < -2 and temperature range of 1.8 < logT [K] < 4.6, while previous 
tables were usually restricted to T > 10 3 K. We find that, while the presence of Deuterium 
does not significantly alter the opacity values, the presence of Lithium gives rise to major 
modifications of the opacities, at some points increasing it by approximately 2 orders of mag- 
nitude relative to pure Hydrogen/Helium opacities. 
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1 INTRODUCTION 

According to the Standard Big Bang Nucleosynthesis (SBBN), the 
most abundant elements in the aftermath of the Big Bang were 
Hydrogen and Helium with very small traces of Deuterium and 
Lithium. These elements make up the so called POP m or primor- 
dial chemical composition. All the heavier elements present today 
had to be produced in the course of the evolution of several gener- 
ations of stars. 

In order to understand the physics of the early Universe, there 
is a need of having appropriate and accurate material functions at 
hand, in particular opacities for a SBBN chemical compositions. 

There is a large body of literature available with opacity cal- 
culations for metal-free (Z=0) H/He mixtures, starting with the 
Paczynski I- VI mixtures (Cox & Tabor 1976, hereafter CT76), the 
calculations of Stahler et al. (1986, hereafter SPS86) and Lenzuni 
et al. (1991, hereafter LCS91), and more recently Harris et al. 
(2004, hereafter HLMT04). 

These calculations were restricted to temperatures between 
10 3 K and 10 4 K. Specific metal-free opacity sets have been cal- 
culated by the OP project (Seaton et al. 1994, 10 3 75 . . . 10 8 K) 
and OPAL (Alexander & Ferguson 1994; Iglesias & Rogers 1996, 
10 3 . . . 10 87 K). More recently, OPAL calculations have been ex- 
tended (Eldridge & Tout 2004) to 10 3 . . . 10 10 K and densities 1 
-8 < log R < 7 with the main application for CNO enhanced 
opacities in stellar evolution. 

All these opacity calculations, however, only considered pure 



1 R = p/T^ is a parameter which replaces the density in order to keep opac- 
ity tables in rectangular format when spanning many decades in temperature 
(for more details see Appendix D); T 6 = T/IQ 6 K 



Hydrogen/Helium mixtures with no metals (mass fraction Z = 0), 
but different ratios X/ Y of the mass fractions of Hydrogen (X) and 
Helium (Y). It has been argued that due to the assumed very small 
abundances of Deuterium and Lithium these element do not play 
any significant role in the opacities and thus in the evolution of 
POP in objects. 

Moreover, while there are good low-temperature opacities 
available for a solar chemical composition (e.g. Semenov et al. 
2003), zero-metallicity opacity tabulations have been restricted to 
temperatures above 1000 K so far. 

In this paper we present calculations of both, Rosseland and 
Planck mean opacities for primordial matter including all the three 
elements (Hydrogen, Helium, and Lithium) including Deuterium 
isotope and present their absorption properties. 

In Sect. 2 we discuss the quantitative composition of primor- 
dial matter at a given temperature and density. We then describe 
the relevant absorption processes (Sect. 3) and present our Rosse- 
land and Planck mean opacities (Sect. 4). Here, we also assess the 
influence of Lithium on the opacity for different Lithium contents. 
Before we compare our results to previous tabulations (Sect. 6), 
we give analytic calculations for the molecule formation timescale 
in Sect. 5 in order to estimate the time needed to reach chemical 
equilibrium which is one of the underlying assumptions in our cal- 
culations. Finally, we summarize our conclusions in Sect. 7. 



2 PRIMORDIAL MATTER 

Primordial Matter as created in the SBBN consists of Hydro- 
gen (H), Helium (He), Deuterium (D), and Lithium (Li). Before 
WMAP(Spergel et al. 2003), the observed abundances of these el- 
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Table 1. SBBN concordance abundances (Coc et al. 2004) from first-year 
WMAP results 

ements have been used to constrain the baryon-to-photon ratio. In 
the framework of a ACDM cosmology WMAP provided the value 
this parameter with high accuracy. While the abundances of H and 
He were fairly well known for some time, the uncertainties in the 
determination of the abundances of D and Li were significantly re- 
duced. Now, the D abundance as derived from cosmological param- 
eters is consistent with direct observations, while the observed Li 
abundance still lacks a factor of 3 compared to the SBBN results. 
This may indicate either systematic effects in the observations or 
new physics (for an in-depth discussion of this issue, see Coc et al. 
2004). We summarize the resultant abundances according to SBBN 
and WMAP in Table 1 (after Coc et al. 2004) and take them as our 
fiducial POP in mixture. For later comparison purposes we define 
also a Li-free mixture in Table 1. 



2.1 Equation of state (EOS) 

We consider 20 species of primordial matter: electrons e~, Hydro- 
gen species H~, H, H + , H 2 , H 2 , and H3, Helium species He, He + , 
He ++ , and He 2 , Deuterium species D, D + , HD, and H 2 D + , and 
Lithium species Li, Li + , Li ++ , LiH and LiH + . 

The computations assume chemical and thermodynamical 
equilibrium. Equilibrium constants (see Table 2) are used when- 
ever they are available in the literature, otherwise appropriate Saha 
equations are taken. 

For the partition functions the most recent H3 partition func- 
tion 2h+ of Neale & Tennyson (1995), for Hj Standi (1994a) 
is used. For HD and H 2 a partition function for the equilibrium 
constant is calculated using the energy levels for HD (Stancil & 
Dalgarno 1997) and 2 H 2 . All other partition functions needed for 
the Saha equations are approximated by the ground state statistical 
weight. 

Combining all the equilibrium constants, there are 5 equations 
utilizing mass conservation of all 4 elements and charge neutrality. 
These equations are solved using an iterative procedure (see Ap- 
pendix A). 

2.2 Equilibrium HJ + H 2 <-> H+ + H 

Since the early work of Patch (1968) it has been clear that Hj is im- 
portant for temperatures of 2000 - 5000 K (Tennyson & Sutcliffe 
1984; Chandra et al. 1991; Sidhu et al. 1992; Neale & Tennyson 
1995) as in this range it is the most abundant positive molecular 

2 Data from http : //www . physast . uga . edu/ugamop/ 



ion (cf. LCS91). The improvement of the understanding of H3 en- 
hanced the reliability and enlarged the high-temperature partition 
function which leads to an increase in abundance. 

Although the maximum abundance never reaches that of H 
and Hj, it influences the equilibrium abundances due to the charge 
neutrality condition. A fit to the latest equilibrium constant using 
the latest partition functions can be found in Appendix C. 

2.3 The influence of D and Li on the EOS 

Deuterium does not affect the equilibrium abundances of the Hy- 
drogen and Helium species given its small abundance and isotopic 
relation to Hydrogen. 

Lithium belongs to the group of Alkali metals. It therefore can 
easily be ionized at comparably low temperatures. Fig. 1 shows 
Lithium providing free electrons at comparably low temperatures 
(T* 800 .. . 3000 K). These enhance the H" but diminish the H+ 
abundance in the high-density limit, while there is no change in the 
H and H3 abundances for lower densities, where only the electron 
abundance is increased. Although being the least abundant species 
in the mixture, it influences the equilibrium abundances of the Hy- 
drogen species. 

Depending on the densities, both elements can form hydrides 
at low temperatures (T ^ 1000 K and T g 300 K, respectively) 



3 ABSORPTION PROCESSES 

In our calculations we take into account 

• Thomson scattering; 

• Rayleigh scattering; 

• Free-free absorption; 

• Bound-free absorption including 

- photoionisation and 

- photodissociation; 

• Bound-bound absorption; and 

• Collision-induced absorption. 

For a detailed list of the processes considered including references, 
see Table 3. 

3.1 Collision-induced absorption (CIA) 

CIA is basically the van der Waals interaction between a pair of 
not neccessarily polar molecules or atoms. The interaction induces 
a dipole momentum which leads to absorption. Although called 
collision-induced absorption, this should not be confused with the 
dipole moments created by an external field. CIA induced dipole 
momenta are in this sense "permanent". The usual timescale of 
this interaction is of the order of nano-seconds, therefore we ex- 
pect broad and diffuse absorption spectra. 

A general theory of CIA has already been developed in the 
1950s (mainly by Van Kranendonk 1957a,b; Van Kranendonk & 
Kiss 1959). The first application to astrophysics has been done by 
Linsky (1969) in the case of late-type stars. During the last decade, 
this theory has been improved by semi-analytical quantum mechan- 
ical computations obtained from first principles using newly avail- 
able dipole moments and an improved line-shape theory (e.g. Bo- 
rysow et al. 1989) 

The importance of CIA has been shown in the model atmo- 
spheres of cool, low-metallicity stars (Borysow et al. 1997) and 
cool white dwarfs (J0rgensen et al. 2000). It is also important in 
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Figure 1. The influence of the presence or absence of Li on the EOS for two densities: (a) p = 10 g cm ~ 3 , POP m; (b) p = \Q- U gem" 3 , Li-free; (c) 
p = 10~ 5 g cm~ 3 , POP m; (d) p = 10 -5 g cm -3 , Li-free. The different lines have the following meanings: solid red line - electron abundance; dotted blue lines 
- abundance of neutral molecules; dashed purble lines - positive ions; dash-dotted black lines - negative ions; dash-triple-dotted green lines - abundance of 
neutral atoms. 



the atmospheres of planets. Observational confirmation has been 
put forward by Bergeron & Leggett (2002) in the case of two ultra- 
cool white dwarfs. It also has been shown to be important in the 
zero-metallicity opacity calculation by Lenzuni et al. (1991) and 
Harris et al. (2004), the first one still partially relying on the Linsky 
spectra. 

In the cool, low-metallicity gas the importance of CIA arises 
because there is no dust absorption. The main absorber is H 2 ab- 
sorbing itself only due to quadrupole transitions which usually are 
very weak compared to dipole transitions. Therefore CIA intro- 
duces the possibility of more powerful absorption in these environ- 
ments. Due to the strong dependence of the absorption coefficient 
on the number density of the contributing species (oc n a n B , see Sect. 
3.4) we expect this absorption to be important in high density re- 
gions. 



3.2 Bound-bound absorption 

Bound-bound absorptions resulting from 8 species have been in- 
cluded in our opacity calculations (see Table 3). While the Einstein 
coefficients and energy levels are available in the literature, the ex- 
act modeling of the lineshape itself proved to be difficult. We tried 



a simple Doppler profile and a Voigt profile (following Shippony & 
Read 1993) including Doppler broadening and radiative damping 
for the Lorenz wings. In both cases the contribution of the molecu- 
lar lines to the Planck opacities execeeds that of the continuum by 
several orders of magnitudes while it does not considerably influ- 
ence the Rosseland mean. 

However, for a more realistic modeling of the lineshapes Stark 
broadening for atomic lines (already considered for atomic H) and 
collisional broadening for the molecular lines (already considered 
for H2 through CIA) have to be taken into account. We consider 
Stark broadening for the first 3 series of atomic Hydrogen (Stehle 
& Hutcheon 1999). This is an important opacity source even for 
the Rosseland mean as it is able to fill the "valleys" between the 
absorption edges of bound-free H absorption in the monochromatic 
absorption coefficients. 

All other bound-bound absorptions (except H3 , see below) are 
treated according to the integrated line absorption coefficient (App. 
B) and are only considered for the Planck mean. We verified this 
approximation by running a extensive computation including Voigt 
profiles for the lines in both Rosseland and Planck mean opacity 
calculation, but never reached differences larger than 10%. 

Given the importance of line contributions, the Planck and 
Rosseland averaging only yields approximative estimates of the 



4 M. Mayer & W.J. Duschl 



H2 <-> H + H 
H~ <-> H + e~ 
H <-> H + + e" 

H+ <-» H + H + 

2 

Hj + H 2 <-> Hj + H 


K _ [H][e-] 
^H" LH -j 

r [H + ][e"] 

ah - [H] 

v [H][H+] 

2 l H 2 + l 

[h+][H] 

3 |H+|[H 2 ] 


Saha (t/ = 4.478 eV) 
Saha ({/ = 0.7556 eV) 
Saha(C7 = 13.6 eV) 
Stancil (1994a) 

Do = 2.6507 eV (Huber & Herzberg 1979) 

see App. C (D = 0.7382 eV, Sidhu et al. 1992) 


HD ^ H + D 

D <-> D + + e" 

B+ + D <-» H 2 D + + H 

+ HD <-> H 2 D + + H 2 


„ m _ [H][D] 

- - [° + ]^] 
K o - [D] 

^ [H 2 D+][H] 
* H 3 + .D [ H3+ j [D] 

_ [H 2 D + ][H 2 ] 
,HD | H +j [HD] 


Saha (U = 4.5167 eV, Stancil & Dalgarno 1997) 
Kd = K H (ass.) 

Sidhu et al. (1992), D = 0.0439 eV 
Sidhu et al. (1992), D = 0.0120 eV 


He «-» He + + e~ 
He + <-> He ++ + e - 

HeJ <-» He + He+ 


^He - [He] 

_ [He ++ ][e-] 
*He+ " [He+] 

[He][Hc+] 


Saha (U = 24.587 eV) 
Saha (U = 54.418 eV) 
Stancil (1994a) 

Do = 2.365 eV (Huber & Herzberg 1979) 


Li «-» Li + + e~ 
Li+ «-» Li ++ + e - 
LiH <-» Li + H 
LiH «-» LiH + + e - 
LiH + <-> Li + + H 
LiH + <-» Li + H+ 


fi Li [Li] 

*Li+ [Li+J 

K .„ , _ [Li][H] 
LiH, 1 - [LiH ] 

[LiH*][e-] 
ALiH,2 - [Lffl] 

^ [Li + ][H] 
*LiH + ,l [Lffl+J 
^ [Li][H+] 
K UH*,2 [LiH + J 


Ku = ^LiH+.l ' ^LiH.l 
Saha (U = 75.640 eV) 
Stancil (1996) 
Stancil (1996) 
Stancil (1996) 
Stancil (1996) 



Table 2. Equilibrium constants 



true monochromatic absorption. Even the line positions are sub- 
ject to change with respect to gas motions, etc. Hence we give two 
sets of Planck opacities (see Section 4). 



3.3 H3 bound-bound absorption 

H3 influences the abundances of the other species in a temperature 
range of 2000 - 5000 K (see Sect. 2.2) as then it is the most abun- 
dant positive molecular ion. 

The bound-bound H3 line absorption has been shown by 
HLMT04 to have an effect of low-metallicity, very low-mass stellar 
evolution. We use the line list of more than 3 million lines (Neale 
et al. 1996), but bin them to frequency intervals of 10 cirr 1 . Our 
tests with different bin sizes ranging from 200 ... 0.1 cirr 1 did not 
produce significantly different results for the Rosseland mean so we 
used this smoothed coefficients which we tabulated in the (v, T)- 
plane. Without loosing accuracy, we save a considerable amount of 
CPU time. Due to the huge number of lines, H3 line absorption ap- 
pears as continuuous absorption. For line absorption in the Planck 
mean, we use the integrated absorption coefficient (see App. B). 



3.4 Database 

For all tabulated absorption coefficients bicubic splines are used un- 
less otherwise stated in Table 3. The tabulated free-free absorption 
were interpolated by bicubic splines, too, but extrapolated to low 
frequencies (cf. John 1994, k vM oc v~ 3 (1 - exp (-hv/k B T)) oc y~ 2 ). 



The CIA data is taken from tables available on the Web 3 ex- 
cept for H 2 /He, where programs available from the same source are 
used to calculate the roto-vibrational and roto-translational spectra. 
These routines are used for temperatures lower than 1000 K, while 
for the high-temperature region (1000 . . . 7000 K) the latest ab ini- 
tio data were applied. Extrapolation to non-tabulated low frequen- 
cies was done according k vCIA oc y 2 . This extrapolation is justified 
as for given temperature at low frequencies, the absorption coeffi- 
cient for CIA(A/B) (e.g. Meyer et al. 1989) 

a{(o) = —n A n B cog((o, T)\l-e 2 *<^ j 

is only proportional to lim^o <o(l-e 2 "'b' ) oc a> while the mainly 
used BC spectral shape (Birnbaum & Cohen 1976; Borysow et al. 
1985) is constant. In the formula above, u> = 2nv is the angular 
frequency, n A and n B are the densities of the species in amagar 1 4 , 
while T denotes the temperature, c the speed of light and h Planck's 
constant. 

3.5 Neglected processes 

There are more CIA data available in the literature (e.g. Borysow 
& Meyer 1988; Gustafsson & Frommhold 2001b, for CIA of 
HD/(He,Ar,H 2 and H)), but their maximum in absorption roughly 

3 see http: //www. astro. ku.dk/~aborysow/programs/ 

4 1 amagat corresponds to a particle number density of 2.686755 ■ 
10 19 cm" 3 
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1000 K < T < 7000 K, 25 cm" 1 < v < 200088 cm" 1 
J0rgensen et al. (2000) 




CIA(H/He) 


H + He + Av -> H + He 


1500 K < T < 10000 K, 50 cm" 1 < v < 1 1000 cm" 1 
Gustafsson & Frommhold (2001a) 




CIA(H 2 /H) 


H 2 + H + Av -> H 2 + H 


1000 K < T < 2500 K, 100 cm" 1 < v < 10000 cm" 1 
Gustafsson & Frommhold (2003) 



Table 3. Scattering and absorption processes 
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Figure 2. Absorption coefficient k s (v) + K a (v) for different densities and temperatures: (a) p = 10 15 g cm 3 , T = 10 2 6 K; (b) p = 10 4 g cm ~ 3 , T = 10 3 K; 
(c) p = 10~ 14 g cm" 3 , T = 10 3 - 2 K; (d) p = 10 -9 g cm" 3 , T = 10 3 6 K. Also shown are the contributions of the different absorption processes. 



coincides with the H 2 /H 2 and H 2 /He CIA in position and height. 
Given the small number fraction of D and the kind of absorption 
(dipole radiation), for the abundances considered here these absorp- 
tion never dominates. 



3.6 Absorption coefficients 

We calculate monochromatic continuum absorption and scattering 
coefficients K a j(v, nj,p, T) and K s j(v,rij,p, T) for all the 28 contin- 
uous absorption and scattering processes (see Table 3) using the 
cross sections cr a/s i (v, T) weighted with the abundance nj(p, T) of 
the contributing species. 

For the absorption processes we correct for stimulated emis- 
sion (factor 1 - exp(-hv/k B T), Mihalas 1978) unless not already 
included in the tabulated absorption coefiicents. We sum them both 
up and divide by the density p to get the mass absorption and scat- 
tering coefiicents K s (v,p, T) and K a (v,p, T) (cf. LCS91) 

^ i 

+ X - Yj <r-ufo r K<A 7>*(A T)(l - (1) 



K,(y,p,T) = -Y^(r sXi {v,T)nj(p,T) 

P i 

+ °W v > T)nj(p, T )"k(p, T) (2) 
^ i 

The indices "1" and "2" indicate 1- and 2-body absorptions. 
Figs. 2 and 3 show sample absorption coefficients. 



4 ROSSELAND AND PLANCK MEAN VALUES 

Using them together with the integrated line absorption coefficients 
(Bl) calculations of the Rosseland and Planck mean opacity are 
done according to 

i = j*_ r 1 dB Vdy 

K R (p, T) 4oT 3 Jo k s + K a dT 
1 C°° \ 1 

K P (p, T)=—fi J o KaBydv + £ Kp M (4) 

where the factors in front of the integrals are the normalization con- 
ditions. 
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Figure 3. Absorption coefficient K a (v) for different densities and temperatures including line absorption: (a) p = 1CT 15 g cm 4 , T = 10 22 K; (b) p = 
lfT 8 g cm" 3 , T = 10 3 - 25 K; (c) p = 10" 6 g cm" 3 , T = 10 4 K; (d) p = 1CT 14 g cm" 3 , T = 10 3 3 K. Also shown are the contributions of the different absorption 
processes. 



The typical range of our frequency grid is 10 . . . 10 7 cm" 1 
which corresponds to photon energies of 1.24-(10" 3 . . . 10 3 ) eV. We 
use typically 20000 logartithmically equidistant grid points with 
additional points added for a better resolution of bound-free ab- 
sorption edges. 

In this contribution, calculations of Planck means are carried 
out for two sets: one set only including the continuum contribution 
and one including all contributions added via the integrated absorp- 
tion coefficient. The only exceptions from this scheme are , CIA 
and Stark broadening of atomic H. The latter obviates opacity cal- 
culation for p > 10" 5 g cm" 3 and T > 10 4 K. 

The two Planck sets are two extreme cases, where the com- 
plete Planck opacity set (continuum+lines) can be used for "zeroth- 
order" calculations, while the Planck continuum is suitable for 
more detailed modeling (including lines depending on the motion 
of the gas, etc.). 



4.1 Rosseland mean 

The upper part of Fig. 4 shows the Rosseland means for our POP 
in Li content f u = 4.15 • 10" 10 and the Li-free composition. 

The p-T-plane can be divided into 4 regions split by the lines 
p ~ 10" 9 g cm" 3 and T ~ 3000 K. At low density scattering 
dominates while true absorption is dominating the high-density re- 



gion. For scattering we distinguish between Rayleigh scattering 
and Thomson scattering (see Figs. 2a+c) in the low- and high- 
temperature regimes. True absorption is divided into mainly bound- 
free and CIA absorption (see Figs. 2d+b). 



In the mid-temperature region there is a steep transition be- 
tween small and large opacity values at low and high temperatures, 
respectively. For the low-density region the steep gradient is simply 
produced by the sharp decline/appearance of H 2 , H and e" (evident 
in the dominance of scattering involving these species) while at 
high densities bound-free absorption of H" is the dominant mecha- 
nism which itself exhibits a strong temperature dependence. 



A comparison of the Li-free with the POP in case shows the 
importance of taking an even minute Li contents into account (as 
was already indicated by the number densities in Fig. 1). The early 
ionisation of Li increases the H" abundance which at high densities 
shows up as an extension of the H" dominated opacity area and as 
an increase in the opacity value. At lower densities the free elec- 
trons increase the opacity due to Thomson scattering. The presence 
of Li, however, decreases the Rosseland means at intermediate den- 
sities due to the destruction of ions, (cf., Figs. 4a and b). In Fig. 
5a we show these differences quantitatively. 
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POP III Lithium-free 




log p [g cm ] log p [g cm ] 



Figure 4. Rosseland (a+b) and Planck mean (without lines: c+d, with lines: e+f) values for the POP III case (fu = 4.15 ■ lO" 10 , left: a,c,e) and the zero-metal 
case (fu = 0, right: b,d,f) 



4.2 Planck mean 



4.2.1 Continuum absorption 



In panels c-f of Fig. 4 we show the Planck mean opacities for the 
Pop m chemical composition (panels c and e) and for the Li-free 
mixture (d and f). In panels c and d we show the values for con- 
tinuum absorption only, while in panels e and f lines are taken into 
account. 



For continuum absorption the Planck means show a rather sim- 
ple three-part structure approximately following the change from 
molecular to atomic to ionic H, respectively. In the molecular H 
region the only dominant true absorption is CIA of H2/H2 pairs, 
whereas for atomic H bound-free absorption of IT and for ionised 
H bound-free absorption of H are the dominant absorption mecha- 
nisms. 

Li influences the opacity in a small stripe at the boundary be- 
tween the CIA (H 2 /H 2 ) and the bound-free IL absorption where 
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the influence of CIA (H/He) and bound-bound of Hj decreases. 
The overall difference is similar to the Rosseland means, although 
the positive deviations are smaller, while the negative ones are 
much larger. At this transition there is a local minimum of the 
Planck mean, as H 2 dissociates/forms, while H and H has not fully 
been formed/destroyed yet. 



4.2.2 Line Absorption 

Taking into account line absorption changes the result consider- 
ably: At regions where H is mostly in its atomic or ionic form, 
bound-bound transitions of atomic H yield the largest contribution 
to the opacity (see Fig. 3c). At lower temperatures, the opacity con- 
tribution of CIA (H 2 /H 2 ) only remains dominant for high densities 
(p > 1CT 8 g cm -3 ), whereas at lower densities line absorption is the 
dominant mechanism (see Figs. 3a,b,d). With line absorption the 
local minimum of the opacity prominently shows up: The drop in 
opacity reaches 7 orders of magnitudes within a very small temper- 
ature interval (at T ~ 2000 K) and low densities. 

D is an important absorber only for temperatures below ss 
100 K, as then the lowest-lying rotational transition of HD at ss 
128 K outperforms the H 2 transition at 510 K in its contribution to 
the Planck mean. 

The direct influence of Li on the Planck means comes into play 
through the absorption of atomic Li (cf. Fig. 3b) and LiH (Li hy- 
dride, see Fig. 3a). Although the number density of these species is 
much smaller than those of D (itself being of considerably smaller 
abundance than H and He), they increase the opacity by one and a 
half orders of magnitude. The indirect influence of the FT absorp- 
tion and H3 destruction is still present. 

Atomic Li absorbs mainly through the 6708 A transition. It 
is stronger than the quadrupole absorption of H 2 through its large 
Einstein coefficient despite taking place at a comparatively long 
wavelength. However the influence of this absorption ceases with 
increasing temperature where Li gets ionized. 

For LiH, the Einstein coefficient is larger than that of H 2 tran- 
sitions due to the presence of a strong dipole moment, but it does 
not reach the values of the 6708 A transition of atomic Li. LiH 
becomes dominant due to the smaller energy differences spacings 
between its rotational levels. This considerably reduces the wave- 
length of the base transition. Thus the Planck function for temper- 
atures lower than 510 K (lowest rotational state of H 2 ) still catches 
rotational and vibrational transitions of LiH (lowest rotational state 
at 21 K) while the H 2 contribution decreases. 

The double-peaked CIA dominated opacities at temperatures 
below 1000 K can be attributed to the the isotropic overlap induced 
dipole and the quadrupolar field induced dipole via isotropic polar- 
izability components in the fundamental band of H 2 /H 2 (Birnbaum 
et al. 1996; Meyer et al. 1989, cf., our Fig. 2). 



4.3 The quantitative influence of Lithium on the opacity 

Here, we summarize the differences the presence or absence of Li 
in the chemical composition of the matter makes. In the Rosseland 
and Planck means Li manifests itself by 

• its early ionisation; 

• the destruction of HJ in its presence; 

• the 6708A bound-bound transition of atomic Li; and by 

• the molecular absorption of LiH. 




-15 -10 -5 

log p [g cm 3 ] 



Figure 5. Difference between POP III (/ Ll = 4.15 ■ 1CT 10 ) and the Lithium- 
free case (fu = 0) for Rosseland and Planck mean opacities: (a) Rosseland 
mean; (b) Planck mean for continuum case; (c) Planck mean with lines 

The first two mechanisms show up in both the Rosseland and 
Planck means, while the latter only appears in the Planck line opac- 
ity. All mechanisms except the HJ destruction lead to an increase 
of the opacity. 

In Fig. 6 we show for different Li fractions (log/ Li = 
-10.3 • • • - 8.3) the differences from Li-free composition. 

The positive deviations (A log k > 0) show a steady increase 
with increasing Li content. The negative deviations (predominantly 
due to destruction of H3 ) increase for Li fractions f Li < 5 • 10~ 10 , but 
are being drowned by the increase of opacity for larger Li fractions. 

5 A MORE DETAILED ASSESSMENT OF THE 
CHEMICAL EQUILIBRIUM 

As we have seen, at low temperatures H 2 is the main absorber. Due 
to the lack of dust which would act as a catalyst for H 2 forma- 
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Figure 6. Deviations from Li-free (fij = 0) opacity for different Li contents. 
A (log k) > indicates positive deviations, while A (log a:) < indicates 
negative deviations relative to Li-free opacity. The fiducial POP in content 
is given including the uncertainties. 



tion (for a comparison, see Glover 2003), H 2 formation and there- 
fore cooling and absorption is small. The main H2 formation chan- 
nels are via H and at low and high temperatures, respectively 
(e.g. Galli & Palla 1998; Lepp et al. 2002). These 2-body reac- 
tions operate on rather long timescales. However at higher densi- 
ties (n > 10 s cnT 3 ) 3-body reactions become efficient (Palla et al. 
1983). 

In the following, we look into the timescales for molecule for- 
mation in order to assess where these chemical timescales exceed 
the hydrodynamical one (for which we take the free-fall time scale). 
Thus we can estimate at which densities and temperatures the as- 
sumption of chemical equilibrium is satisfied at all. 

For all timescale calculations we assume that initially the con- 
stituents of the molecules are mainly in their neutral atomic form 
and subsequently get converted. This approach can be justified 
given the fact that in the cosmological fluid H and He are in their 
neutral atomic form starting from a redshift of z = 1 100 while for 
Li we take z ~ 100 (Galli & Palla 1998; Lepp et al. 2002). 



5.1 Molecular Hydrogen H 2 

The H 2 formation timescale r H , due to the reaction H + H + H 
H 2 + H can be estimated using k 4 of (Palla et al. 1983) to 

-2 



t h , = 2420yr-( ) 

2 y \300K/ 



(5) 



, 10" 14 g cm 

In contrast to the rapid 3-body-reaction the only relevant ionic 




log p [g cm" J ] 

Figure 7. Comparison of free-fall timescale with the molecule formation 
timescales at T = 300 K and a H? fraction of 10 -5 



H 2 formation channel at low temperatures is via H while at higher 
temperatures H 2 is formed via H 2 (McDowell 1961; Peebles & 
Dicke 1968; Saslaw & Zipoy 1967). The H~-consuming reaction 
is H + H — * H 2 + e~. This reaction depends on the H abundance 
which never is the most abundant H species and rapidly reaches 
chemical equilibrium as H abundance determining reactions pro- 
ceed on timescales much shorter than those for the H 2 chemistry 
(Abel et al. 1997). Thus taking into account the H 2 formation rate 
for the ionic reaction channel is not advisable, as it would only give 
the timescale on which all H is depleted. The maximum H 2 abun- 
dance due to this reaction never exceeds the available H density. 

Following arguments similarly shown by Bromm & Loeb 
(2004), there exists a critical density (as a function of temperature), 
at which H 2 formation proceeds within one free-fall time 



Pcrit.H, 



3.2 • irr 



: 425 yr • 



10- 



gcm 



V300K/ 



g cm 



5.2 Deuterium hydride HD 



HD is formed mainly via deuteron exchange with H 2 (H 2 + D — > 
HD + H) For the timescale we use rate coefficient (1) of Galli & 
Palla (2002) and get for the timescale (setting T = 300 K) 



^hd = 1.97 • 10 3 yr • 



10- 14 g cm- 3 



fu 2 
10- 5 



Note that thd depends on the H 2 fraction. Thus with increasing H 2 
abundance this timescale is reduced. 



5.3 Lithium hydride LiH 

For LiH we use the low-temperature rate of reaction (20) of Stancil 
et al. (1996) and get 

"I . rr, . -0.28 



Tuh = 140 yr . 

1 10 14 g cm- 3 



(300 k) 



and 



, j .-0.56 

^uh = 1.08.10 -«(_) g, 
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The situation, however, will change, once there are data available 
of LiH producing 3-body reactions (see discussion in Stancil et al. 
1996). This will modify the density dependence of the formation 
timescale (t oc p~ 2 ). The same applies for HD. 

At densities p > 10~' 4 g cnr 3 chemical equilibrium of H 2 is 
reached within one free-fall time scale or less. Thus, for our the 
computational domain (1(T 16 g cnr 3 < p < 1CT 2 g cnr 3 ) the as- 
sumption of chemical equilibrium is justified, with the exception 
of the two lowest-density decades (which we tabulated neverthe- 
less for computational convenience). LiH can be formed even more 
rapidly and thus can be regarded as being in equilibrium for densi- 
ties p > 1CT 15 g cnr 3 . HD is rather hard to form, as long as there is 
no significant H 2 abundance because the formation timescale cru- 
cially depends on the H 2 fraction. However, we have seen in Sect. 
4 LiH to be a significant absorber while HD is not. Thus the HD 
abundance is of only minor relevance for the opacities presented. 

t H d, together with the other timescales discussed in this Sec- 
tion, is shown in Fig. 7. 



6 COMPARISON WITH EXISTING CALCULATIONS 

In order verify our results, we compare the Rosseland means of our 
calculation with existing calculations of Z = opacities by setting 
fu = /d = 0. We restrict the comparison to more recent opacities 
(Section 6.1-6.4). A common difference of all these opacities com- 
pared to our data can be found for temperatures around 10 3 2 " 3 4 K. 
This comes about due to our use of the newly calculated Hj equi- 
librium constant in conjunction with the symmetry factor in the H 2 
dissociation equilibrium. 

An overview of previous low metallicity opacities is shown in 
Table 4. 



6.1 Stahler et al. (1986) 

The match with the SPS86 opacities (Fig. 8a) is good as long as 
H~ is the dominant absorption mechanism. For low temperatures, 
the deviations are due to the use of CIA data of Linsky (1969) and 
Patch (1971) while we use the more recent Borysow data. The devi- 
ations at higher temperatures and densities are due to the influence 
of the Stark broadening of H lines which have been neglected in the 
SPS86 calculation. 



6.2 Lenzuni et al. (1991) 

The comparison with Lenzuni et al. (1991) is done for both Rosse- 
land and Planck means (Fig. 8b&c). 

At higher temperatures and densities the deviations in the 
Rosseland mean are again due to Stark broadened H lines. For 
lower temperatures the deviations are only in the 20-30 % range, 
much smaller than with respect to SPS86. This is due to their use 
of the partially newly available CIA data of H 2 /H 2 and H 2 /He 
collisions for the roto- vibrational and roto-translational transitions 
while still having to use extrapolations for the overtones Linsky 
(using data by 1969). 

For the Planck means the differences are considerably higher. 
This is due to the linearity of Planck averaging which prefers peaks 
in the opacity (CIA at low temperatures, H lines at higher tempera- 
tures and lower densities). 



6.3 Alexander & Ferguson (1994) & Iglesias & Rogers (1996) 
Z=0 sets 

The comparison to the OPAL opacities for the Grevesse & Noel 
chemical composition at Z=0 shows a much better concordance 
over a larger temperature range (Fig. 8d&e). 

In the high-temperature limit (logT > 3.8) the Rosseland 
mean opacities are consistent at the 10% level (Stark broadening). 
For the very highest temperatures, deviations are due to the approx- 
imation of the statistical weight of H by the ground state which no 
longer is true. At lower temperatures we again have a 10-20% de- 
viation due to CIA. 

While the inclusion of Stark broadening for the Rosseland 
means shows a good consistence with the OPAL data, there are 
deviations in the Planck means. These differences are due to our 
coarse frequency grid which does not correctly trace all the peaks 
of the Hydrogen lines. 

6.4 Harris et al. (2004) 

The comparison to the HLMT04 data in Fig. 8f shows an overall 
good agreement with our data. The discussion is analogous to the 
OPAL comparison, except that now we are consistent at the 10% 
level accuracy for lower temperatures due to their use of H/He, 
H 2 /H and more recent H 2 /H 2 and H 2 /He CIA data. 

At densities p > 10~ 7 g cnr 3 and temperatures log T > 

3.5 there are, however, important differences. These are due to 
HLMT04's use of two different approaches to calculate the 
abundance in order to overcome machine accuracy problems. In 
their H 2 dominated regime, they take a dissociation energy of 
4.52 eV whereas we take 4.478 eV (following Huber & Herzberg 
1979). In our calculations, we do not find any problems regarding 
machine accuracy. The second reason for the deviations is probably 
due to the difference in H 2 absorption coefficients where we use the 
Stancil (1994b) data while they use Lebedev et al. (2000). 



7 CONCLUSIONS 

Rosseland and Planck mean opacities for primordial matter have 
been calculated using the most recent available data for the absorp- 
tion mechanisms. 

We tabulate the opacity data in Tables El, E2 and E3, 
for a temperature range 1.8 < log(T/K) < 4.6 and -16 < 
log(p/g cirT 3 ) < -2 for our POP in matter composition (cf. Ta- 
ble 1). These are the first POP m opacity tables for temperatures 
T < 1000 K. 

In order to make application as easy as possible, we provide 
two sets of Planck opacities: Planck means including only contin- 
uum absorption and those including molecular lines. Higher reso- 
lution tables including routines for bicubic interpolation tables are 
available from the authors upon request. 

It has been shown (see Sect. 4 and 4.3) that the small number 
fraction of Li leads to a significant change in the opacity values at 
temperatures T < 4000 K compared to a pure H/He mixture. The 
differences can reach up to 2 orders of magnitudes. There are four 
processes which change the opacity: 

As an alkali metal, Li gets ionized at comparatively low tem- 
peratures. It therefore increases the number of available electrons 
considerably. They increase the Thomson scattering contribution 
to the Rosseland mean at low densities while increasing the H~ 
absorption at higher densities. For both the Planck continuum and 
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Figure 8. Comparison of our Rosseland and Planck mean opacities with previously published Z = 0-opacities: (a) kr by Stahler et al. (1986) (X = 0.72, 
¥ = 0.28, and Z = 0); (b) & (c) kr and a:r by LCS91 (X = 0.70, ¥ = 0.30); (d) & (e) kr and kj> by Alexander & Ferguson (1994) for log T < 4 (log T < 4.1 
for Planck mean) and Iglesias & Rogers (1996) for log T > 4 (X = 0.70, Y = 0.30); (f) a:r by Harris et al. (2004) (X = 0.72, Y = 0.28). The contour lines 
correspond to differences of the opacity in %. Dashed lines indicate negative, solid lines positive differences relative to our opacities. The increment between 
neighbouring contours is 10 %. 
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Reference 


Opacity type 


Density range 


Temperature range 


Chemical composition 


Cox & Tabor (1976) Paczynski I-IV mixture 


KR 


-12 < logp < -4 


1500 < T < 12000 K 


X = 0.0... 1.0 


Stahler et al. (1986) 


kr 


-15 < logp < -3 


1500 < T < 7000 K 


X = 0.72 


Lenzuni et al. (1991) 


kr/p 


-12 < logp < -0.3 


1000 < T < 7000 K 


X = 0.72 


Alexander & Ferguson (1994) 


kr/p 


-7 < log/? < 1 


3.0 < log 7" < 4.1 


X = 0.0... 0.8 


Seaton et al. (1994) OP Project 


kr/p 


-7 < logR < -1 


3.5 < log T < 8.0 


X = 0.0 . . . 0.7875 


Iglesias & Rogers (1996) OPAL 


kr 


-8 < logfi < 1 


3.75 < log 7/ < 8.70 


X = 0.0... 1.0 


Harris et al. (2004) 


kr 


-14 < logp < -2 


1000 < T < 9000 K 


X = 0.0... 1.0 


This work 


kr/p 


-16 < logp < -2 


1.8 < log 7/ < 4.6 


X = 0.75 



Table 4. Characteristics of Z = opacity calculations or opacity calculations with a Z = subset (p in g cm 3 ;R = pT 6 3 ; T in K; T 6 = T/W 6 K) 



line case only the H absorption at low densities (no scattering) is 
enhanced. 

In the presence of metals (e.g. Lithium) HJ is being destroyed. 
As Hj is the most abundant positive ion at temperatures around 
3000 K and densities p > 10" 12 g cnr 3 it influences the opacity 
indirectly in changing the chemical equilibrium when being de- 
stroyed. Furthermore, in part it shows up directly via bound-bound 
transitions of HJ at intermediate densities. For higher densities this 
is being drowned by the indirect increase of the H~-absorpfion. 

The Planck means including line absorption are changed due 
to absorption by atomic Li and molecular LiH. 

For atomic Li the most important transition is the 6708A fea- 
ture. The Planck function, however, is most sensible to features at 
6708A for temperatures of approx. 7000 K, a temperature at which 
Li is ionized, at the latest. Absorption via this transition is still im- 
portant at temperatures around 1500 K because the corresponding 
Einstein coefficient is very large compared to the quadrupole tran- 
sitions of molecular and atomic H. 

At temperatures T < 500 K the influence of H 2 is ceasing as 
the lowest lying transition of H 2 has an equivalent temperature of 
512 K. For lower temperatures either HD or LiH contribute to the 
opacities. LiH, however, has got a larger dipole moment and thus 
has got much larger Einstein coefficients and transitions at much 
lower temperatures (as low as 15 K). 

A critical point in primordial chemistry is the formation of 
molecules. Formation times have been calculated for H 2 , HD and 
LiH. For densities p > 10~ 14 g cnr 3 molecule formation proceeds 
within one free-fall time (except HD, which does not play any sig- 
nificant role in the POP m case). Hence the calculation is valid for 
densities larger than that, provided the free-fall time scale being the 
shortest timescale relevant, other than the chemical. We give values 
for densities as low as 10~ 16 g cm" 3 for numerical convenience. 

In comparison to previous calculations we find a good agree- 
ment of our results when neglecting the newly added absorption 
mechanisms and the contributions of Li. 

Based on our new opacities the influence of Li on the different 
stages of POP in star formation and evolution can now be assessed. 
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APPENDIX A: EOS FOR POP III 

Primordial matter is assumed to consist of a number fraction /e of 
element E. Depending on the density, we define the number densi- 
ties of the involved elements 

,, _ (/h, 1 - /h - /d - /li, /d, fu] p 

■W|H,H(;,D,Li) - -. 77 77 — -7 

4 - 3/ H - 2/ D + f u m p 

We need 4 equations (A1-A4) to describe the conservation of the 
number densities of the elements (mass conservation) and an addi- 
tional one (A5) to take care of charge neutrality. The number den- 
sity of species X is called [X]. 



+ [H] + [H + ] + 3[H 3 + ] + 2 ([H 2 ] + [H 2 + ]) = 


= N H 


(Al) 


[He] + [He + ] + [He ++ ] + 2[He+] = 


~- N Hs 


(A2) 


[HD] + [D] + [D + ] + [H 2 D + ] = 


~- N D 


(A3) 


[LiH] + [Li] + [Li + ] + [Li ++ ] = 


~- N u 


(A4) 



[e-] + [H-]-[H + ]-[H 2 + ]-[Hj]-[He + ] 

-2[He 2 + ] - [D + ] - [H 2 D + ] - [Li + ] - [Li ++ ] = (A5) 



Note that we neglect the H contribution of the D and L species in 
the H sum. Given the relative number fraction of D and Li with 
respect to H this assumption is satisfied. The equilibrium constants 
are given in Table 2. 



Rosseland and Planck mean opacities for primordial matter 1 5 



For the A «-> B+C Saha equilibria we use equations of the 



form 



n K(T) ■■ 



with £ on = 2nm;k B /h 2 = 1.80 • 10 10 K cm- 2 and f dlss = 
2mn p kB/h 2 = 3.30 • 10 13 KT 1 cm -2 . U is the characteristic energy 
for ionisation or dissociation, Q the partition function of the con- 
tributing species. For H 2 -dissociation we use £ diss and correct the 
equilibrium constant with the symmetry factor for a homonuclear 
molecule (Tatum 1966). 

For the solution of the system of equations we define a critical 
Temperature r crit = 200 K above which we solve the full system in- 
cluding the charge neutrality condition. We call this the ionic limit. 
Below r cril we neglect all ions and solve (A1-A4) in the molecular 
limit. 

In the ionic limit, we start with an initial guess of [e~], solve 
eqns. (A1-A4) to derive a new [e~] from the neutrality condition 
eq. (A5). Iteration is done with [e~] — > y[e~] new [e~] oM . After 20- 
50 iterations convergence is reached. The solution of eqns. (A1-A4) 
is done analytically and would involve a cubic equation in the case 
of eq. (Al), otherwise only quadratic equations are involved. The 
cubic term in the H species equation arises due the inclusion of the 
abundance. As this species never is the main contributor to the 
number density, we neglect it in the number density conservation. 
Note that we also neglect the H containing D and Li species in the 
H conservation equation. Since we have /d//h <k 1 and fulfil <K 1 
this assumption is justified. 

In the molecular limit we have to solve analytically only eqns. 
(A1-A4) neglecting ions. 



APPENDIX B: THE INTEGRATED ABSORPTION 
COEFFICIENT 

The absorption coefficient a v [ cirT 1 ] for the 2 — > 1 transition can 
be written as 



c 2 / -51l\ 

a v = ^j*(v)niA 12 11 - e ^ I , 



where O(v) the profile function, >i[ is the number density of atoms 
in state 1, A l2 the Einstein coefficient and v\ 2 the transition fre- 
quency. In relation to the ground state we have 

B\ C 2 I - £ 1 - E 2 \ 

a y = 2i — O(v)n A , 2 <f *5T _ e " ¥T , 
#o 8ttv 2 V / 

where E\ and E 2 are the energies of the states relative to the 
ground state. In terms of the total number density n = = 
(no/go) Zi ft exp(-£, A B D = n Q(T)/g 

a v = 8l 8^km* {v)nAl2 ( e ~^ ~ e ~^)- 

Via giA 12 = g 2 A 21 we easily can transform to the usually used 
Einstein coefficient for spontaneous emission, A 21 . 

Evaluating avdv yields (J[ 0(v)<iv = 1) the integrated ab- 
sorption coefficient Xv a [ crrr ' s ~'] 



Xvi2 ~ 8i 



%nv 2 Q(T) 



nA l2 \e W - e W\ 



where we made use of the approximation v ~ const across the line. 
For Planck averaging we furthermore assume B v ss const 



across the line. Then the line contribution of the 2 — > 1 transition is 

1 f y a v B v dv 



P f v B Y dv 



1 



I a v B, 



crpT* • 

"V 



n ^ 



1 



crpT\^ 6 ' Snk 2 Q(T) 



A 12 le W - e ^ T )S V 



hen 



4no-pT 4 Q(T) 



ftcyj_,2 

- 1 



(Bl) 



where we introduced the wavenumber k^ 2 corresponding to the 
1 — > 2 transition. Multiplied with p/n, the last expression only de- 
pends on the temperature. Calculated once for each species sum- 
ming up all the lines for the temperature grid chosen, it can be 
added directly to the Planck mean with the density weight n/p. 

Note that we make use of the Einstein coefficient A [2 , which 
is related to the Einstein coefficient for spontaneous emission A 2l 
through the Einstein relation giA i2 = g 2 A 2l . 



APPENDIX C: EQUILIBRIUM + H 2 <-> + H 

We have fitted the equilibrium constant log, K to the 6 parameter 
function (x = log, (77 K)) 

5 



io gl0 *(r) = 



with the coefficients 



7079 K 



1 + Y aj 



a x = +(0.365793 ±0.02124) 
a 2 = - (0.624287 ± 0.03488) 
a 3 = + (0.388540 ± 0.02095) 
a 4 = -(0.102057 ±0.00546) 
a 5 = + (0.940462 ± 0.05203) • 10~ 2 

This fit causes maximum deviations of 4 % of A" for 300 K < T < 
7000 K, whereas the maximum deviation never exceeds 8 % below 
10 4 K. 7079 K correspond to the enthalpy of 16300 K for the reac- 
tion HJ + H 2 — » Hj + H, divided by the log 10 factor due to basis 
conversion. 



APPENDIX D: A REMARK ON THE STANDARD 
FORMAT OF OPACITY TABLES 

Opacity mean values are usually tabulated in the (R, T)-plane, 
where R = A.T 6 is the temperature in units of 10 6 K, p the density. 

T 6 

Using this ^-parameter, the tables remain rectangular over large 
temperature ranges. 

We have the gas (P g ) and radiation (P r ) pressure (here in the 
optically thick case) 



P g =P 



kJ_ 

pm p 



3c 



(Dl) 



Calculating the ratio and taking the logarithm we have 
^ = l0g (S) = 1 ° 8 (4^) + 3f -^ 22 - 5 + 3f ^ (D2) 
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Subtracting the logarithmized definition of R = p + 18 - 3T, 
we have 

f- if = 4.5 

For R = -4.5 gas and radiation pressure are equal. As the 
range of R usually tabulated is -8 < R < 1 (OPAL, OP Project), 
R = -4.5 is in the middle of this coordinate range. 

The physical meaning of introducing R is that one remains in 
the physical regime where gas and radiation pressure compete with 
each other. Too large departures from this equality require changes 
in the EOS. 

Therefore the each column of opacity tables in the style of 
OPAL & OP represent different gas/radiation pressure ratios. 



APPENDIX E: SAMPLE TABLES 

Sample tables for Rosseland and Planck means are given in Ta- 
bles El, E2 and E3 (will be available in the accepted version). 



